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expressed in relation to four dimensional perturbation theory are in exact agreement with the 
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are determined to six loops and 0(l/Nj). The computation of the exponents in the Schwinger 
Dyson approach is also provided and an expression for rj in arbitrary covariant gauge is given. 
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1 Introduction. 



Recent developments in the area of multiloop calculations in quantum chromodynamics, (QCD), 
have included the determination of the MS four loop f3- function, [Q] , and the four loop quark mass 
anomalous dimension, [Q, 3- These calculations involved current state of the art analytic and 
symbolic manipulation techniques in the evaluation of the order of 10, 000 Feynman diagrams. 
Such high order calculation are necessary to refine our understanding, for example, of the running 
of the quark masses, [||, g|. Indeed such four loop calculations have built on the early lower order 
work of several decades of |Q, ||, ||, 0, g] for the QCD /3-function and M, [U^, 11] for the quark 
mass anomalous dimension. Whilst we believe the results of Q and [^, [3[] will not be superseded 
by a five loop calculation for quite some time, it is possible to probe the higher order structure 
of the QCD renormalization group, (RG), functions using techniques alternative to conventional 
perturbation theory. One such approach is the large Nf expansion where Nf is the number of 
quark flavours. In this method the Feynman diagrams contributing to the determination of an 
RG function are reordered according to the number of quark loops and evaluated by considering 
at leading order those graphs which are simply a chain of bubbles. The next order is represented 
by those bubble chains which have either one vertex or self energy correction and so on. Clearly 
this large Nf reordering will cover the information already contained in explicit perturbative 
loop calculations but will reveal new information beyond the current orders. For simple scalar 
field theories such as the O(N) c/> 4 theory or the O(N) nonlinear a model, the large N technique 



has been developed to determine information on the RG functions at 0(1/A 3 ), |12| , 13]. That 
impressive programme of Vasil'ev et al exploited ideas from critical renormalization group theory 
and determined rather than the RG functions themselves, the related critical exponents. These 
correspond to the RG functions evaluated at the d-dimensional fixed point of the /3-function 
which in O(N) c/> 4 theory is the Wilson-Fisher fixed point. The techniques developed in the early 
work on scalar theories, [12, 14], have been applied to QCD at leading order in 1/Nf in [15, pl|. 
Further, information on the anomalous dimension of the twist-2 operators fundamental to the 
operator product expansion used in deep inelastic scattering has been produced as a function 
of the operator moment, [^|. Although it is important that these calculations are extended to 
next order in 1/Nf, it turns out that there are several issues which need to be addressed. First, 
no rigorous proof exists of the critical equivalence of QCD and the so called nonabelian Thirring 
model, (NATM), which underlies all the large Nf computations. The original observation of 
the connection between the models in |l8j was effectively at 0(1/Nf). Second, the analytic 
regularization commonly used in 1/N calculations breaks the gauge invariance of the theory. 
As far as we are aware all other regularizations spoil the masslessness of the propagators and 
this therefore makes higher order calculations virtually impossible. We note that dimensional 
regularization, which is ordinarily used in conventional perturbation theory, is not applicable in 
large Nf work since the theory remains divergent in arbitrary dimensions. This second obstacle 
is much more serious and might seem to destroy the possibility of developing a sensible 1/Nf 
scheme. The resolution of the issues of the proof of the critical equivalence of QCD and the 
NATM as well as demonstrating that the critical exponents evaluated in the latter model using 
non-gauge invariant regularizations do in fact match those of QCD are of extreme importance 
and represent the central results of this paper. By way of application and in order to support 
the scheme which will be developed in practical calculations we examine the quark and mass 
anomalous dimensions and determine them at new order in the large expansion which is 
0(1/Nj). A preliminary version of our results was given in [19]. 

There are various motivations for examining these two RG functions. First, as in conven- 
tional perturbation theory one always needs to deduce the wave function renormalization of the 
fundamental fields of a theory before considering the renormalization of the other parameters 



2 



and operators in a Lagrangian. Likewise in the large Nf approach the evaluation of the wave 
function exponent, 77, needs to be determined first. Indeed in the context of the phenomenolog- 
ical application of the large Nf technique the programme of evaluating anomalous dimensions 
of the twist- 2 operators of deep inelastic scattering at 0(1/Nj) cannot proceed prior to the 
determination of 77 at the same order. Second, although rj is fundamental it depends on the 
covariant gauge parameter introduced in the gauge fixing of the NATM Lagrangian. Conse- 
quently, it is not as fundamental a quantity as a gauge independent exponent. Therefore, to 
further understand the large Nf method in QCD at 0(1/Nj) we have undertaken to consider 
the gauge independent quark mass anomalous dimension. Indeed as part of our calculations we 
address the issue of the choice of gauge in the large Nf method. Whilst it may seem that the 
evaluation of two exponents at 0(1/Nj) is a formidable task, it turns out that the computation 
of the quark mass dimension does not represent a significant amount of additional calculation 
as there is, in fact, a close relation between the fermionic graphs for the wave function and mass 
dimension which we will establish. 



It is worth noting that our detailed calculation relies on the novel method of |2C, 21] for 



the computation of critical exponents in 1/N expansion. In comparison with the standard 



technique [12| based on the direct solution of the Schwinger-Dyson (SD) equations the latter 
approach allows us to express the critical exponents through the corresponding renormalization 
constants whose explicit form is very similar to those computed using dimensional regulariza- 
tion. The immediate advantage is that it allows one to use techniques similar to the infra-red 
rearrangement of conventional high order perturbative calculations to evaluate the three loop 
critical point Feynman diagrams needed for the present calculation. This approach, which was 
originally developed to analyse the 0(N) nonlinear a model appears to be universal, and can 
be easily extended to the case we will consider. Together with some calculational shortcuts it 
makes the evaluation of almost all the relevant graphs quite straightforward and simplifies com- 
putations significantly. In addition we will also present the Schwinger-Dyson formalism for the 
determination of r] at 0(1/ N?) in QCD. This builds naturally on the earlier large Nf calcula- 



tions in quantum electrodynamics, (QED), |22] , 23, 24]. In those papers the original technique of 



[12] in solving the SD equations was followed and in extending the formalism to QCD here will 
provide an interested reader with a unique opportunity to compare both approaches. Finally, 
it is worth stressing that the consistency of our final results for the quark and mass anomalous 
dimensions at 0(1/Nj) with the corresponding four loop perturbative results of || 1C, [Tl|, ||, ||] 
provides a nontrivial check on the validity of either approach. 

The paper is organised as follows. Section 2 is devoted to introducing background formalism 
and providing the proof of the critical equivalence of QCD and NATM. In section 3 we develop the 
large Nf method explicitly in the context of QCD and the NATM. The details of the evaluation 
of the three loop diagrams to determine the quark anomalous dimension are described at length 
in section 4. Section 5 is devoted to the development of the Schwinger Dyson formalism to 
determine the quark anomalous dimension. The computation of the extra three loop graphs 
needed to determine the mass anomalous dimension is given in section 6. The final results for 
the quark anomalous dimension and mass anomalous dimension are collected together in section 
7 where we also derive new information on the coefficients of the respective RG functions. Several 
appendices contain results which were fundamental to our calculations. 
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2 Background. 



The QCD Lagrangian in d = 4 — 2e dimensional Euclidean space reads 

C = ^ipf 1 + ^F« v F a ^ + ^(d-A) 2 + d,c a (D»c) a , (2.1) 

where ip is the quark field belonging to the fundamental representation of the colour group, 
1 < I < Nf, A® is the gluon field, c a and c a are the ghost fields in the adjoint representation of 
the colour group, £ is the covariant gauge parameter and g is the coupling constant. The field 
strength tensor F^ v and the covariant derivative are defined by 

= d,A a v - d v Al + f abc A^Al 

and 

D„ = (d^-iA^) , 

where T a are the colour group generators in the corresponding representation and f abc are the 
structure constants with [T a ,T b ] = if abc T c . To ensure the coupling constant, g, is dimensionless 
below four dimensions we rescale it in the standard way by setting g — > M e g, where the parameter 
M has dimensions of mass. The partition functions of QCD are defined as 

(0 1 (x 1 )...O n (x n )) = Z~ l J D<$> 1 (x 1 )...O n (x n )e W {-S} , (2.2) 

where $ = {A, *0, ip, c, c} is the set of fundamental fields, Oj(xj) represent a basic field or a 
composite operator and Z is a normalizing factor. As usual, the divergences arising in the 
calculation of (|2.2| ) are removed by the renormalization procedure at each order of perturbation 
theory. Namely, provided that the averaging in ( |2.2| ) is carried out with the renormalized action 
all correlators of the fields taken at different space-time points will be finite. The renormalized 
action S l ^($,e), (e = {g, £}), has the form of the action ( |2.1| ) with the fields and parameters 
being replaced by the bare ones 

S R ($,e) = S($ ,e ) , $o = Z$ $ , e = Z e e . 

Here Z$ = {Za, Z^, Z c } and Z e = {Z g ,Z^} are the renormalization constants. However this 
procedure does not guarantee the finiteness of the Green functions with operator insertions, 
which contain additional divergences. To remove these extra divergences one should renormalize 
the composite operators as well. In general, the renormalized operator reads 

[Oi] R = J2 Z *kOk , (2.3) 
k 

where the operators have canonical dimensions equal to or less than those of the original 
operator Oj, and Z^ is the mixing matrix of renormalization constants. An operator is called 
multiplicatively renormalized if the matrix Z^ is diagonal, giving [Oi]n = ZiOi where there is 
no summation over i. 

The renormalization group equation (RGE) for the one-particle irreducible n-point Green 
function with the insertion of k multiplicatively renormalizable composite operators reads 

(m<9 m + p g d g + (3^ - n $7$ + J2 TO^ T(xi, . . . , x n+k ,M, g,£) = (2.4) 
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where we use the shorthand notation n$7$ = n A j A + n^j^ + n c j c . The RGE functions /3 g , 
(3^, 7$, 'yoi which are the respective beta functions and the anomalous dimensions of the fields 
and composite operators, are defined as follows 



09 



M^-g(M) 



7* 



7o, 



d 

dM Ul 



(2.5) 



Formally, the parameters g and £ enter ( |2,4| ) on the same footing and should be considered as 
independent charges. However, the gauge fixing paremeter ^ is introduced into the theory by 
hand and cannot enter an expression for a physical quantity. Thus for proper (gauge invariant) 
objects the term P^d^ drops out of ( |2.4D , which then takes the form of an RGE of a single charge 
theory. 

Our subsequent analysis relies heavily on the existence of a non-trivial infra-red, (IR), stable 
fixed point g* of the ci-dimensional /^-function, P g (g*) = 0, for large values of Nf. The /3-function 
has been calculated in MS using dimensional regularization up to 0(a 6 ) terms, where a = (g/2-7r) 2 
= a s /ir, in Q, 6, [?], |8|. We only record the first few terms here, while the full four loop result 
can be found in 



A* (a) 



(d-4)a + 



11 



\t f n, - l {c A 

79 



a 2 + 



1 



l -C F T F N f + \c A T F N f - X lc\ 



12 



205 



-C F T F Nf + m C A Jpf/ + —C F T F Nf - —C F C A T F N f 



864 

from which it follows that 
3e 1 



432 

1415^ 2m „ 2857 



1728 



a 



16 

1 + 0(a 5 ) 



288 



a* 



T F Nf 



+ -^pq ( 33C< Ae - [27C F + AbC A ] e 2 




(2.6) 



(2.7) 



The Casimirs for a general classical Lie group are defined by 



TV T a T 



T F S ab 



acd fbcd 



c F i , r ca r 



C A 8 



ab 



(2. 



It immediately follows from ( |2.4|) that the Green functions of gauge invariant operators are 
scale invariant at the critical point g*. In other words 



G(Xxi 



\ Di G(xi) 



where Di is the scaling dimension of the corresponding Green function. Moreover, due to the 
IR nature of the fixed point, this index determines the power of the leading term of the IR 
asymptotic behaviour of the Green functions (pi — > 0) near the critical points and this is stable 



against the perturbation of the action by IR irrelevant operators [25|. On the other hand Green 
functions of gauge dependent objects, such as the propagators of the basic fields which will in 
general depend on £, are not invariant under scale transformations. Although one may restrict 
attention from the outset to gauge independent quantities, since they have physical meaning, 
it is also possible and convenient to choose £ so that all Green functions are scale invariant. 
Evidently, this is equivalent to the condition /%(#*,£*) = 0. Since 



2£(e + 7A + Pg/g) 



(2.9) 



one concludes that the equation /%(<?*, £*) = has two solutions. One is £;* = whilst the other 
is 7,4(9*) £*) = — £■ Bearing in mind that our main aim is the development of the 1/iVf expansion 
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we choose the first solution, £ = 0, since the latter gives £ ~ Nf, which leads to problems in 
the construction of the 1/Nf scheme. The origin of the above two solutions for £ becomes more 
transparent if one tries to write down the most general form of the gluon propagator satisfying 
the requirements of both scale and gauge invariance. Indeed, scale invariance yields 



M 



(p 2 ) 



-( 



(2.10) 



where P^ v and P^ v are the transverse and longitudinal projectors, respectively, and A and B 
are constants. As is well known, p5| , radiative corrections do not contribute to the longitudinal 
part of gluon propagator. Hence, (?JL = t\g 2 M 2e p\ v p~ 2 . This implies, that if a 1 then £ must 
vanish, £ = 0. On the other hand for £ ^ then one must have a = 1 which is easy to check 
is equivalent to ■ja = — e corresponding to the canonical dimension of the field. Earlier work 
concerning the relation of scaling and conformal symmetry in the context of gauge theories has 
been given in [26, 27]. We also note that in an abelian gauge theory the condition 7^4 = — e 
follows directly from the Ward identities that implies the scale invariance of all Green functions 
at the critical point g = g* in any gauge [24]. 



As is well known from the theory of the critical phenomena [£5[ physical systems which look 
quite different at the microscopic level may exhibit the same behaviour at the phase transition 
point. An example of this universality is the fixed point relation between the Heisenberg ferro- 
magnet and (i-dimensional (ft 4 field theory. In what follows we construct the theory belonging to 
the same universality class as large QCD but which has a simpler structure. To this end we 
develop the 1/Nf expansion correlators of the type given in ( |2.2| ) and analyse their behaviour 
in the IR region. As usual, the first step is to integrate over the fermion fields in the functional 
integral to obtain the effective action for the gluon field 



Sf EE Nf 



trlnO-i/ a T a ) + 



M 



-2c 



4g 2 



+ 



M 



-2c 



-(OA) 2 + d^cf (D»c) a 



(2.11) 



where bearing in mind that g 2 ~ 1/Nf we have set g 2 = g 2 /Nf. Assuming that Ay- is a large 
parameter, one can evaluate the integral with action ( 2.11| ) by the saddle point method that 
generates the systematic expansion for the correlators. 



We can now examine which terms in (|2.11| ) contribute to the leading IR asymptotics of the 
correlators. We start our analysis with the gluon propagator in the Landau gauge, £ = 0, which 
in leading order in 1/Nf is 



GW(p) 



n; 



2\d/2-l 



Af 



-2c 



-1 



P, 




P, 



(2.12) 

Since 2e = 4 — d > the contribution to the propagator from the second term in the brackets 
of ( pull) , {F/i v ) 2 , is suppressed by (p/M) 2e in the IR region as compared with the contribution 
form the fermion loop which is the first term in the brackets. Thus the asymptotic form of the 
gluon propagator in the IR limit is fully determined by the fermion loop contribution 



G^{p) 



IR 



LIU 



(2.13) 



Careful analysis shows that diagrams with triple and quartic gluon vertices do not contribute to 
the leading IR asymptotics of the correlators either. To see this we note at first that the term 
(F^) 2 in the action ( ^TT|) may be considered as an ultraviolet regulator with M playing the role 
of a cutoff. Indeed, due to this term the gluon propagator decreases as 1/p 2 as p — > 00, which 
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makes all diagrams convergent. We note that the appearence of divergences of the type A 2 ^ is 
prohibited by gauge invariance. Next, if we let G denote an arbitrary diagram, then rescaling 
the variables in the momentum integrals, l{ — > Ml' if one finds that up to some prefactor M fce , 
the parameter M enters the integrand together with the external momenta only. In other words 
G(q ex t,M) = M ke G{q ex t/M). Therefore the limit q ext — ► corresponds to the limit M — > oo, 
which can be regarded as the removal of the regularization. Further, since the most IR singular 
contributions we are interested in come from the integration over the region where l\ ~ 0, one 
can replace the gluon propagator by its asymptotic form fl2.13|) . Then it can be easily checked 
by simple dimensional analysis that the leading IR singularities come from the diagrams without 
gluon self-interaction vertices. On the formal level this follows from the fact that term (F^ u ) 2 
is accompanied by the factor M~ 2e and vanishes in the M — » oo limit. So one concludes that 
this term does not influence the critical properties of the theory and according to the general 
scheme |25|] should be excluded from the action. Therefore we obtain the theory defined by the 
action 

S nat M = 4>(0-ifl a T a )4; + ^(rr^dAf + a„*W + f abc d^<fAic c , (2.14) 

It, 

which in the Landau gauge has the same critical behaviour as QCD to all orders in the 1/iVf 
expansion. Of course, for gauge independent quantities it is true in any gauge. We have adapted 
the usual form of the gauge fixing condition to ensure that the transverse and longitudinal parts 
of the gluon propagator have the same momentum dependence. It is interesting to note that in 
order to arrive at ( [2.14 ) one can start from the theory with manifest gauge invariance which is 



determined by the action S = ip($ — i^L a T a ) , ip with ghost and gauge fixing terms in turn arising 
from the application of the Faddeev-Popov procedure to the functional integral. 



Power counting shows that the theory ( 2.14 ) is renormalizable within the 1/Nt expansion 
and the renormalized action has the form 

S R = ^ - iZ 2 ^fi a T a ^ + ^{U- e / 2 dAf + Zzd^d^c + Z A f abc d^cA b ^c c (2.15) 

where we assume, of course, that a gauge invariant regularization is used. Due to the Slavnov- 
Taylor identities the renormalization constants Z{ are related by 

Z 2 Z^ = Z A Z^ . (2.16) 

This equation was used in the exponent formulation to determine the ghost anomalous dimension 
at 0(1/Nf) in the Landau gauge in [jT5|j . As was proved above in the Landau gauge, the critical 
properties of QCD and this new theory which we shall refer to as the non-abelian Thirring model, 
(NATM), are identical. Therefore one can use the NATM to deduce the QCD RG functions. 
This equivalence at leading order in 1/Nf was noted in |l8j and used to deduce various exponents 
at 0(1/ N f ) in g [2§, g|, The extension of these calculations to 0(1 /Nf) requires 

special care. The main one is the necessity of using a gauge invariant regularization which was not 
crucial at 0(1/NA. The conventional dimensional regularization is not applicable here, since 
the gluon propagator behaves as (p 2 ) 1_d / 2 and the theory remains logarithmically divergent 
in any dimension d. To our knowledge most other invariant regularizations such as higher 
derivatives spoil the masslessness of the propagators, which makes higher order calculations 



virtually impossible. Usually in 1/N calculations the analytical regularization of [12] is used. 
However, this breaks gauge invariance. 

We now consider how to reconcile gauge invariance with the calculational advantage of using 
massless propagators. First, we break the gauge invariance of ( 2. 14| ) from the beginning by 



introducing a new coupling, A, for the ghost-gluon vertex of (2.14). For A = 1 we recover the 
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original model but the theory remains renormalizable for arbitrary A as well. The only effect 
will be that the identity ( |2.16| ) will no longer hold. The bare coupling Ao is related to the 
renormalized one, A, by 

A = Z\\ = Z 4 Z 1 Z 2 1 Z^ 1 X, (2.17) 



where the Zi will now also depend on A. Let us suppose that we used an invariant method, 
such as regularization by higher derivatives [29], to regularize this extended theory. Then it 
immediately follows from ( 2.16| ) and ( 2.17 ) that the equality A = 1 for the renormalized coupling 
implies that Ao = 1 as well. Therefore we conclude that A = 1 is a fixed point, /?a(1) = 0. The 
existence of this fixed point is the key point and it does not depend on the regularization used. 
So using any other regularization can only change the position of the fixed point with in general 
A* = 1 + 0(1/Nf). What is important, though, is that the anomalous dimensions calculated 
at the critical point, 7(A*), are scheme independent and, hence, coincide with the anomalous 



dimensions deduced in the original model (2.14). Therefore one can use the regularization which 
is most convenient from the computational point of view. Moreover, since we do not need to 
consider diagrams with external ghost legs, then the only diagrams depending on A are those 
with a ghost loop. As is evident from counting powers of 1/Nf these are themselves 0(1/N^). 
So at this order it is sufficient to set A = 1. 



3 Methods. 

As was shown above the RG functions of QCD at the critical points can be deduced from 
the more simple NATM. The use of the A-regularization 12, |TJ] allows us to retain massless 



propagators which means that the calculation of the higher order Feynman integrals can be 
achieved. However, the price to be paid for this is the loss of the property of multiplicative 
renormalizability. Indeed, formally, the regularized action has the form 

S A = $frp - i$4 a T a i() + 2 ^ 2A (l3~ (e ~ A)/2 0M°) 2 + <9 M c a ^c a + \f abc d^c a Alc c 
- l -Al{x)M-^K^(x,y)Al(y) + -A^x)^ (x,y)A a v (y) . (3.1) 
We use K^ v {x,y) = K IJ,p {x — y) for the inverse fermion loop in coordinate space with 

5^(^-2?). M 

where n = Nf Tp Tr sp j nor I. The regularized kernel K^(x,y) is defined by 

K%{x) = C{A)K^(x)\x\ 2A (3.3) 



and the constant C(A) = 1 + 0(A) will be specified below. The last term in (|3.1|) forms part 
of the interaction, while the penultimate one is assigned to the free part of the action. Since the 



last two terms in the action (3.1) are non-local, they are not renormalized and therefore the full 



renormalized action takes form 

si = z$w> - iz 2 iprT a iP + (□- (e - A)/2 aM«) 2 + z 3 d^c a 

+ Z A Xf abc d^Ay - \Al{x)M-^K%{x,y)Al{y) + ^(x)^(z, y)A a v (y) . 



(3.4) 



S 



Obviously, Sj± cannot be brought into the form ( |3.lD by the redefinition of the fields and con- 
stants. Consequently, the theory with action (|3.1|) is not multiplicatively renormalizable. A 
more detailed discussion of this topic can be found in the [14, 30]. The absence of multiplicative 



renormalizability means the following. We recall that in multiplicatively renormalizable theories 
the Green functions calculated within two different subtraction schemes are related to each other 
as follows 

G n \%i i ■ ■ • j x n > 9) = A n G n (x\ , . . . , x n , Zg g) , 

where the coefficient A n = ZQ and Z$ is the field renormalization constant. In the case under 
consideration this equality holds as well. However, A n ^ Z 1 ^ any longer with A n depending on 
n in a nontrivial way. This of course does not contradict the statement about the IR equiv- 
alence with QCD, since the IR asymptotics remains unchanged. Nevertheless, the absence of 
multiplicative renormalizability prevents us from using the standard method of RG analysis for 
the determination of critical exponents. The more widely used method for this is the method of 
self-consistency equations, which is based on the direct solution of the Schwinger-Dyson (SD) 
equations, (l^, |l4|, Pqj . Another approach has been developed recently in [20, 21| and will 
be briefly discussed below in the context of QCD. 



3.1 Extended model. 



The convenience of the RG method from the computational point of view consists of the possi- 



bility of calculating critical exponents via the renormalizaton constants given in (2J3). To retain 
this calculational advantage, following the lines of [14, ffifl , we restore the multiplicative 
renormalizability of the model in question by attaching two additional charges, u and v, to the 
last two terms of the action (|3.1[). Namely, we consider the model with the action 



u 



Suv = S NATM - ^(x)M' 2A K^(x,y)A a M + -^(^(^^(y) . (3.5) 
Obviously, the initial model ( |3,1| ) is recovered by the special choice of the parameters u = v = 1. 



(1 




(1 




Figure 1: The effective gluon propagator for u, v 7^ 1. 

As distinct from the initial model the propagator of the gluon field in the model ( |3.5| ) with u 
and v couplings has a more complicated form. Indeed, at v ^ 1 the last term in the action ( |3.5| ) 
does not ensure the exact cancellation of the simple fermion loop insertion in the gluon line and 
one should sum up all such insertions (see figure |]), that yield for the transverse part of the 
gluon propagator 



G ± {p;u,v) 



-gHp) ( 1 + { — 

u \ u 



(3.6) 



with t = C(A)(M 2 /p 2 ). This theory is obviously renormalized. In addition to the standard 
renormalization constants, Z\, . . ., Z4, which depend now on the couplings u and v, one should 
add two new ones Z u and Z v to take account of the renormalization of the couplings u and v, 



u 



uM- 2A Z„ 



vZ v , 



Zv 



Z 



-2 
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Being multiplicatively renormalizable this model can be analysed with standard RG methods. 
In particular the basic RG functions are 

7$ = Md M In Z$ , f3 u = Md M u = 2u (A - j A ) , (3 V = Md M v = - 2v~/ A , 
Pz = Md M t = - 2$ (A - 7 a) , Px = M8 M \ , (3.7) 

and 

Z^ = z\^ 2 , Z c = z\j 2 , Za = Z2 Z^ 1 , Z\ = Zi Z z 1 Z\ Z 2 1 , Z% = Z 2 A . (3.8) 

Before discussing the question of the relation of the RG functions of the extended model to the 
critical exponents of QCD, we derive a compact expression for the RG functions at the point 
u = v = 1. Henceforth we will adopt the minimal subtraction (MS) scheme in which the general 
renormalization constants have the form 



00 

z = 1 + £^ (n W,£,A) . 



n=0 A 

Then taking into account the finiteness of the RG functions and using ( |3.7|) one obtains 

7 = M8 M In Z\ u=v=1 = 2(ud u -^)Z^(u,v,^\)\ u=v=1 . (3.9) 

Since there are no derivatives with to respect A and v, we can set A = A* = 1 + 0(1/Nf) and 
v = 1 from the very beginning. In this case only the first term in ( |3.6| ) survives and, therefore, 
the operation {ud u — £c^) applied to a diagram gives simply the number of the gluon lines, ua, 
in the latter 

(ud u -^)T(u,v,^X)\u=v=i = - rur(l,l,£,A) . (3.10) 
The formula (^^) can be rewritten as follows 

7 = - 2]T4 r) 4 1) (l,U,A) , (3.11) 
r 

where the sum runs over the set of all diagrams. Thus to calculate the RG functions defined by 
one can put u = v = 1 from the very beginning but take into account the number of gluon 



lines in each individual graph. To represent ( 3.11 ) in a compact form we attach the factor g to 
the gluon propagator 

G„ v gG^ , (3.12) 

so that each diagram with k internal gluon lines acquires a factor g k . Then for the anomalous 
dimensions of the basic fields, ( j3~S| ) takes the form 

7$ = - 2gd g Z^\ g=i , (3.13) 

where the renormlization constants Z$ are calculated within the model with u = v = 1 and the 
modified gluon propagator ( 3.12j ). 

Similarly, the matrix of anomalous dimensions of a system of composite operators, {Oj}, 
which mix under renormalization, is given by 

7i fc = 2gd g z[l ] + 5 ik n k ^<y§ , (3.14) 

where the mixing matrix, Zi k , is defined in the standard way from the condition for the Green 
functions of the renormalized operators O^, 

Of = Z ik O k , (3.15) 
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to be finite. Again, Z$ is the coefficient of the simple pole in A. One can easily note the obvious 
resemblance of the formulas (3.13) and (3.14) with those used in the MS scheme in dimensional 
regularization. Therefore using this approach gives a simple and transparent way to compute 
the critical exponents of operators and, importantly, all the machinery of e expansions can be 
adapted for the case in question as well. 



3.2 Critical exponents and RG functions. 



In the previous subsection we have developed the effective tool for the computation of the RG 
functions ^fig = 1rg{ u i v)\ u=v= \ in the NATM. The non-trivial exercise is to show that these RG 
functions, coincide with the corresponding critical exponents 7 cr jj, which determine the scaling 
properties of the correlators. Generally speaking, this is not so and ^rg la-it- Indeed, the RG 
equation for 1-irreducible n-point Green functions with an operator insertion in the extended 
uv model reads 

(\Md M + Pudu + + (3\d\ - n* 7 $] 6 ik + f£ G )r k ( Pl , . . . , Pn+1 ;u, v,£,\) = . (3.16) 

Bearing in mind the critical equivalence between QCD and the NATM in the Landau gauge, we 
set £ = and A = A* and do not display their explicit dependence in the Green functions. Since 
= 0x(^*) = from (|3.7D, the above equation takes the form 



([Md M + Pudu ~ n$7*] $ ih + iRG^kiPU ■ ■ ■ ,Pn+r,u, A) = . (3.17) 



Evidently, due to the presence of the term Pudjj m ( 3.17 ) the latter does not describe the scaling 
properties of the Green function Tk(pi, . . . ,p n+ i;u,v) under scale transformation even at the 
point u = v = 1. To make this more transparent, taking into account (|3.7j ), we rewrite ( |3.17| ) 
as follows 



([ 



Md M ~ n*7*] $ ik + lR G Fk({pi}) 



25 ik 7A (d u + d v )r k ({ Pl }; U ,v) 



u=v=l 



(3.18) 



It is evident that the matrix can be considered as the matrix of the anomalous dimensions 
only if the right side of ( 3.1 8| ) is equal to zero. However, since in the Landau gauge 7^ ^ the 
right side does not vanish identically. This means that knowledge of the matrix 7}^. alone is 
not sufficient for the calculation of the critical exponents. 

The problem of the relation of the ^rg and jcrit has been analysed in |p0| , 21 where it had 
been shown for the example of the 0(N) nonlinear sigma model that the difference 7^ — j C rit 
is 0(l/N 3 ). Thus, up to 0(1/N 2 ) the formulae ( 3.13| ) and (|3.14|) give the true answer for the 
critical exponents. To prove this statement one has to show that the right side of ( p. 18 ) is 
0(1/N 3 ). Since the anomalous dimension of the gluon field is 0(1/Nf), it is sufficient to prove 
that the quantity (d u + d v )Tk({pi}) u, v) vanishes at OCl/Nt) for any Green function. The 

. . u=v=l 



proof given in [ 20 1 is reduced to the problem of checking this property and without any change 



can be simply adapted to the present model. More details can be found in [2C, 21]. Thus we 
have justified the use of formulas fl3,13j ) and ( 3.14 ) for the computation of critical exponents in 
the NATM or QCD in the Landau gauge up to 0(1/ Nf). 

It is worth considering what happens in other gauges. Of course, as has been stressed 
earlier, in this case it is only sensible to consider the anomalous dimensions of gauge invariant 
operators. In general, the anomalous dimensions of gauge invariant operators, such as ipif) 
or (F^ u ) 2 , calculated with ( 3.13j ) and ( 3.14 ) in the extended model, will be gauge dependent. 
Indeed, we broke the gauge invariance of the model explicitly by introducing the coupling A 
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for the ghost-gluon vertex and implicitly by the regularization we have introduced. So it is 
not surprising that the correlators of gauge invariant objects in this model appear to be gauge 
dependent. However, it follows immediately from the critical equivalence of QCD and the 
NATM that at the critical point A = A* all dependence on £ in the correlators of gauge invariant 
operators factorizes 

T(xi,0 = A r (Of(x<) , (3.19) 

where T(xi) is already independent of £. Then the RG equation for the correlator of two gauge 
invariant operators 0\(x) and 02{y) can be cast in the form 

(Md M + 70l + 7o 2 + 7*) f (x, y) = , (3.20) 

where 7* = ffcd^ lnA(£). Bearing in mind that the explicit breaking of gauge invariance occurs 
at 0(1/N f ) only, (A* = 1 + 0(1/N f )), one can deduce that 7* is 0{1/Nf). This immediately 
leads us to the conclusion that up to 0(1/Nj) the RG functions of the gauge invariant operator 
given by the formulas ( 3.13j ) and ( |3.14| ) do not depend on the gauge fixing parameter. 



Finally, since (3 U and V as well as /?£ are proportional to 7a from (|3.7|), one can nullify all 
these beta functions simultaneously choosing the value of the gauge fixing parameter from the 
condition 

MU) _ „ , U _ _ fifcil + o (£) . (3.21) 
In this gauge the RG equation in the extended uv model takes the desired form 

[Md M - n* 7$ ] Sue + J% G )r k ({ Pl }) = , (3.22) 



from which it follows that the RG functions ^rg are the genuine anomalous dimensions 7cHt' 
We conclude this section by summarizing the main results: 

• In the Landau gauge the critical exponents of QCD at the critical point can be calculated 
with standard RG methods via the renormalization constants in the extended NATM with 



the help of the formulas ([□]§) and ( |3H| ) up to 0(1/N? 



f 



• Up to 0(1 /Nj?) the correspondence between the critical exponents of gauge invariant op- 
erators in QCD and NATM holds in any gauge. 

• For the special value of the gauge fixing parameter £ = £4 the formulae ( |3.13j ) and ( |3.14| ) 
yield the anomalous dimensions of the gauge invariant operators in all orders of 1/Nf 
expansion. 

In the next sections we will demonstrate the effectiveness of the above approach by calculating 
the anomalous dimensions of the quark field and the mass operator at second order in the 
1/Nf expansion. The explicit values of the renormalization constants Zi and the basic RG 
functions at 0{1/Nf) are given in appendix B. 



4 Calculation of graphs. 

In this secton we discuss the technical details of the calculation of the diagrams relevant for 
the determination of the anomalous quark dimension. The technique for the evaluation of 
massless Feynman diagrams is widely documented. (See, for example, [[H], 0). Nevertheless, 
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(a) (b) (c) 




(d) (e) 

Figure 2: Diagrams contributing to the computation of rj2- The first graph represents the gluon 
self energy diagrams of figure [|. 

for completeness, we provide a list of basic formulae, such as the rules for the integration of chains 
of propagators for scalar, spinor and vector fields and the uniqueness relation for the different 
type of vertices in appendix A. In what follows we shall concentrate mainly on the calculational 
problems specific to the case under consideration. 

First, we consider the diagrams with external quark legs which are illustrated in figure |2[ It 
transpires that these are not as complicated to evaluate as may appear at first sight. Indeed 
diagrams (a) and (b) are trivial in that they are equivalent to simple chain graphs in the language 
of |0] . Although (c) has been considered in the case of QED we re-evaluate it here in the context 
of the method used to determine diagrams (d) and (e) where the latter differ from each other 
by a colour factor. 

First, for (c) and (d) it is convenient to choose the flow of the external momentum p in 
both diagrams to be along the fermion line connecting the external vertices. The diagrams are 
linearly divergent. Differentiating them with respect to the external momentum p^ one obtains 
a set of logarithmically divergent diagrams, which differ from the initial ones by the insertion of 
a new vertex (— 7«) in each fermion line. This is due to the result 

g i> + i = _ i> + i i> + i , 41 -j 

QpV (p + q) 2 (p + q) 2 M (p + q) 2 

We recall that after subtraction of the divergent subgraphs, in other words after the application 
of the 1Z' operation, the residues of the poles in A in the diagrams do not depend on the 
external momenta. Therefore, we are free to choose an arbitrary route for the flow of the 
external momenta in the each of the resulting diagrams to simplify the subsequent calculations. 
Our choice is the following. For diagram (c) with the insertion of the vertex (—7^) in the vertical 
fermion line we direct the external momenta flow from the upper to the bottom vertices. When 
the insertion is in the right or left fermion line we choose the route of the external momenta flow 
from this new vertex to be along the fermion line nearest to the incoming (outgoing) vertex. 
After this rearrangement all diagrams are easily integrated since they are reduced to simple 
chains. For diagram (d) we choose the inserted vertex adjacent to the external one as the new 
external vertex, and direct the momentum along the fermion line connecting them. Then, the 
resulting diagram is also reduced to chains but with the insertion of the two loop master diagram 
in one of the lines, see figure ||. The evaluation of this two loop diagram is straightforward. 

We now turn to the discussion on the calculation of the diagrams contributing to the gluon 
self-energy. These are shown in figure ||| where (b), (c) and the graph involving ghosts, (d), 
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q = 




p p 

Figure 3: New external momenta routing for the diagram Fig. ^j, (d). 



are again trivial to determine and do not deserve further comment. On the contrary graphs 
(a), (e) and (f) are tougher to evaluate and we give details of their determination. As (e) 
and (f) are related by the same up to the colour factor we will focus on the former. Again 
(a) was evaluated in the QED calculation but we reconsider it here due to the fact that 
new techniques were required to evaluate (e) which can more easily be appreciated in a two 
loop topology. We recall first the standard method for calculating the diagrams with these 

(a) (b) (c) (d) 





(e) 





W666666515 
(f) 




Figure 4: The diagrams contributing to the gluon self-energy at 0(1/ 



topologies in the example of the nonlinear a model, 12]. It will allow us to avoid the un- 
necessary complications related to the nontrivial 7-matrix structure of these diagrams in the 
case of QCD. So, we first consider diagram (a) of figure ||, where now all the lines are as- 
sumed to be scalar ones, with respective indices (2 — A) and (fx — 1) for the wavy and arrowed 
lines. The diagram is superficially convergent but has two divergent subgraphs. If one could 
set A = then the diagram could be easily integrated due to the uniqueness of the 3-point 
vertex. To determine this diagram it was suggested in [12] to subtract and add diagrams 
which have the same singularities as the initial one, but which can be explicitly computed. 






Since the singularities cancel in the difference one can take the limit A — > 0. Then all diagrams 
can be calculated due to uniqueness of the 3-point vertex. The extension of this approach to 
the case of QCD is possible but leads to the certain problems. The first one is the increase 
in the number of diagrams to be evaluated. In other words five instead of the one for the 
nonlinear cr-model. The second and more difficult one is that for A = the quark-gluon vertex 
is not unique. Using the terminology of ]j3~ll l it is one step from uniqueness, which suggests the 
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application of the integration by parts method. This also increases the number of integrals to 
be considered. The calculation of diagrams with propagators having nontrivial tensor structure 
is more involved compared to scalar ones. Therefore it is desirable to try and keep the number 
of diagrams during a calculation to a minimum. The approach which is the most economical 
in this sense is likely to be the one advocated in [31|. Again, we explain the idea in the simple 



example of the a model. We shift the indices of the upper lines in diagram (a) by an amount 
e and the lower ones by — e. In other words we set fj, — 1 — > fj, — 1 + e in the top lines. It 
is important to note that here e is a temporary regularization and ought not to be confused 
with the parameter which is conventionally used as the dimensional regularization in standard 
perturbative calculations. The analytic expression for such a diagram can now be written as 

^F(A,e), (4.2) 

where F(A, e) is a regular function in the vicinity of the point A = e = 0. What is important 
is that we have shifted the indices of the lines in such a way that the pole structure of the 
diagram has remained unchanged. Next, we need to know this diagram up to constant term in 
A at e = 0. Due to the obvious symmetry e — > — e the function F only depends on e 2 , and its 
expansion near A = e = is 

F(A, e) = F + Aiq + A 2 F 3 + e 2 F 4 + 0(A 3 ,Ae 2 ,e 4 ) . (4.3) 

Since the first two terms of the expansion we are interested in are independent of e, one is 
allowed to evaluate the diagram at any value of e. It is convenient to set e = A/2 which results 
in the uniqueness of the upper vertex and the diagram is immediately integrable. It is worth 
noting that one has only to calculate a single diagram instead of five. Applying this procedure 
twice to the diagram with topology (e) one reproduces the known answer ]l2| with a minimum 
amount of calculation. For further discussion it will be important that one is able to determine 
the subsequent, A 2 term in the expansion (|4.3| ) as well. Indeed, bearing in mind that the pole 
term .F(0, e)/A is fully determined by the counterterm diagrams or, equivalently, is given by the 
application of the 1Z' operation to the graph in question, one can easily deduce the coefficient 
i<4. This in its turn allows us to fix the value of the coefficient F$. 

It now seems reasonable to apply this idea to the QCD diagrams. However, one obstacle 
remains in that the quark-gluon vertex is not unique for the choice of exponents we are restricted 
to. Though when the propagator of the vector fields has the conformal form 

G- f (x) = j^y " , (4.4) 

then the 3-point fermion vector vertex is unique if the sum of the vertex indices, 2a + j3, is equal 
to d + 1, where a is the index of the fermion line. For this unique vertex the relation given in 
appendix A holds. However the propagator of the gluon field does not have a conformal form. 
Indeed, at first order in the 1/Nf expansion in the Landau gauge it reads 



G 



, 2(1 - A) x^x v 
^ + (2n - 3 + 2A) x 2 



(4.5) 



where G is some constant. Of course, bearing in mind that the calculation of the diagrams with 
a longitudinal gluon propagator is rather trivial one may represent Q4.5|) in the following form 



GV(x) = A(»,A)G™f(x) + B{^A)Gl v (x) , (4.6) 

where A and B are some constants and g\ v (x) is purely longitudinal in momentum space. In 
such a decomposition it is easy to check that the constants A and B are singular as A — > 0, 
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A ~ B ~ 1/A. The reason for this is that when A = the conformal propagator ( |4.4|) co- 
incides with the longitudinal one. If these constants were finite then one could calculate the 
diagram for the conformal and longitudinal part of the propagator separately. The diagrams 
with longitudinal propagator would be trivial to integrate, whilst for those with a conformal 
propagator could be evaluated with the methods we discussed earlier. Now, the singularity in 
the coefficients of A and B lead to additional difficulties. Namely, due to the additional factor 
1/A in the coefficients A and B one should evaluate this diagram with higher accuracy in A. 
However it is a reasonable price to pay for the considerable simplification which arises from the 
uniqueness of the triple vertex. 

In the following we shall mainly discuss diagram (e) of figure |I| since its evaluation was the 
most difficult. We will focus our discussion on the more important points as the intermediate 
steps are not difficult to reproduce. However, to give an impression about the effectiveness of 
the suggested appproach it is instructive to discuss the calculation of the QED type, diagram 
figure |](a), first. Using the decomposition (|4.6|) one reduces the problem to the calculation of 
the diagram with the conformal propagator ( [4.4|) with index 7 = 1 — A, where the calculation 
of the diagram with longitudinal propagator is trivial. Due to the singularity of the coefficient 
A(A) this diagram should be evaluated up to terms linear in A. Following the scheme discussed 
above in the case of the scalar diagram we introduce the additional regularizing parameter ±e 
in the upper and lower lines. Then the analytic expression for this diagram is given by the same 



formula (4.3), where all the necessary Lorentz indices are implied. For e = A/2 the diagram 
can be integrated due to the uniqueness of the upper vertex, while the coefficient F4, as was 
explained earlier, can be extracted from the one loop conterterm diagrams. Thus, we have 
reduced the determination of the photon self energy diagram to the calculation of four chains. 

We now consider the calculation of the three loop gluon self energy diagram. First, since 



1 - A 




Figure 5: Three loop gluon self energy diagram with e and A regularizations. 

the longitudinal part of the gluon polarization tensor U pu (p) can be easily calculated, it is suffi- 
cient to calculate the above diagram with contracted indices, Hf, which simplifies the algebra 
significantly. Moreover, the contributions of all the diagrams in figure ^ to the longitudinal part 
of the gluon polarization cancel identically. Further, since the calculation of the diagram with 
the longitudinal propagator is trivial, we shall discuss only those with the conformal ones. For 
this diagram we again introduce the additional regularization e in the fermion lines as shown in 
figure ||. The diagram, being superficially convergent with two divergent two loop subgraphs, 
the analytic expression for it is given as follows 

n(A,e) = 1 (n + An 1 + A 2 n 2 + A 3 n3 + e 2 n4 + e 2 An 5 + O(A 4 ,A 2 e 2 ,e 4 )) . (4.7) 

We recall that our purpose is to determine the function n(A, 0) up to the A 3 terms. In other 
words we are interested in the coefficients Hq, III, II2 and 
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(a) (b) (c) 

Figure 6: The three resulting diagrams for G con f <8> G con f case. 

The strategy of the calculation is the following. For e = A one can exploit the uniqueness 
of the lower right and upper left vertices to evaluate the diagram, which allows us to find 
the coefficients IIo and IT and the combination of coefficients II2 + II4 and 113 + 115. To 
extract II2 and II3 we evaluate the coefficients II4 and II5. The former is determined by the 
counterterm diagrams and its calculation is straightforward, while the evaluation of the latter 
is more nontrivial and will be discussed below. At first, however, we consider in more detail 
the computation of the diagram for a special value of the parameter e. We take e = A, which 
restores the uniqueness of the two vertices. 

Using the uniqueness relation one can express the initial diagram as the sum of two loop 
diagrams, which after some algebra can be reduced to the ones shown in figure |6| and a chain 
integral. Here the wavy line with labels v and p denotes the conformal propagator G c °™^ . The 
dashed line with label p and index a is used for the propagator x p /(x 2 ) a . In diagram (c) the 
vertical double line is used to represent ^ ® where each ^ enters from a different fermion cycle. 
We discuss briefly the calculation of each of the diagrams shown in figure ||. 

• Diagram (a). 

This diagram enters the expansion with a coefficient proportional to 7 — 1 = —A so we 
need to calculate it up to 0(A). To this end we shift the indices of upper fermion line 
by e, and the lower ones by — e. The analytic structure of this modified diagram is the 
following 

n a (A,e) = ^ + R'(A,e) = I(fl + e 2 i? 1 + O( e 4 )) + R' Q + R[A + 0(A 2 ,e 2 ) , (4.8) 

where we took into account the fact that due to the e — ► — e symmetry the diagram 
depends on e 2 only. Further, the diagram can be integrated at e = A due to uniqueness of 
the Gross-Neveu type at the upper vertex. Again, the coefficient R± is determined by the 
counterterm diagrams and can be easily evaluated. This is sufficient to determine II a (A, 0) 
with the required accuracy. 

• Diagram (b). 

At first we note that this diagram is finite as A — > and we need it up to 0(A 2 ). Again, 
we shift the indices of the wavy lines by +e and — e. Then the expansion of this new 
diagram in A and e reads 

n 6 (A,e) = AR + A 2 fli + e 2 R 2 + ... . (4.9) 
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When A = the diagram can be calculated exactly which allows us to find the coefficient 
i?2. Next, it can be checked that the diagram can be integrated for e = 2 A as well. Indeed, 
representing {r\ vp — 2x u x p /x 2 ){z — x) p as 

z 2 (z — x) 2 

where we choose the coordinates of the left and right vertices to be and z respectively, and 
the upper vertex to be x, one obtains three diagrams. The first two of them can be calcu- 
lated due to the uniqueness of the upper vertex, whilst the last is a simple chain. We note 
that one needs to introduce an additional regularization 5, by for example shifting the index 
of the lower dashed line of the original diagram by 5, to make each of the resulting diagrams 
finite. Of course, the sum of the diagrams is finite in the 5^0 limit. Again it is sufficient 
to determine ILj(A,0) with the required accuracy, ILj(A,0) = ILj(A,2A) — 115(0, 2A). 

• Diagram (c). 

The last diagram enters the expansion with the coefficient (7 — l) 2 = A 2 . So one needs 
to calculate it with 0(1) accuracy. The calculation runs along the same lines as above. 
We shift the indices of the upper lines by +A, and the lower ones by —A, which evidently 
does not influence the A terms. To simplify the 7-matrix structure of the diagrams it is 
convenient to represent the tensor product of the 7-matrix, 7^ (g) / j p where each 7-matrix 
enters from different traces, as ^ p G p \ <8> G xp 7^, where x is the coordinate of the upper 
vertex and G p \ is the numerator of the conformal propagator flj.4|) . We repeat this exercise 
for 7ct ® 7 CT and *y u ® 7". Then, using the following identities 

tYG pX (x) = - 7A^ , G Xa {x)Y = - ilxil^ 2 

both traces can be cast into the form suitable for the inversion transformation. After the 
inversion the upper vertex decouples from the base (left vertex) and the diagram can be 
immediately integrated. 

The evaluation of the coefficient II5 is grounded on the following observation. Since in the 
momentum space representation the conformal propagator takes the form 

G^W - J^TTE (*V> ~ 2(M _t + A) P± w) < 4 ' 10 > 

the contributions to the e 2 A term come only from the momentum integral involving two longitu- 
dinal propagators pH ®P", or longitudinal and transverse ones, P" (g'P - '". Moreover, because the 
transverse propagator enters the expansion ( |4.10| ) with an additional factor A, the contribution 
from the latter can be extracted again from the consideration of the counterterm diagrams alone, 
and does not require much work. Thus the nontrivial part of the calculation is the determina- 
tion of the diagram in question with the longitudinal propagators. We shall use the momentum 
representation and at the first step contract one momentum from the numerator of each longi- 
tudinal gluon propagators into the fermion traces. After this the initial diagram reduces to four 
integrals. Two of these are identical whilst another is a simple chain integral. Therefore, we are 
left with two nontrivial diagrams which are shown in figure |7[ For later convenience we replace 
the regulator e in the right fermionic triangles by 5. We consider the left hand diagram first and 
denote it by G(A, e, <5). We look for the e 2 term in the expansion of function G(A,e,e). It is 
easy to see that for e — > with A and 5 fixed that G — > as well. Indeed in the e — > limit the 
left hand triangle turns into a fermion loop, which is transverse, so that when it is contracted 
with the incoming momenta it vanishes. Thus one concludes that G(A, e, 5) = e5G(A, e, 5) and 
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H + A 




Figure 7: Two diagrams contributing to the evaluation of H5. The indices of lines are given in 
the momentum space version. 

for our purpose it is therefore sufficient to calculate e5G(0, 0, 0), where we take into account the 
finiteness of the diagram. This implies the following method for its evaluation. We ascribe the 
label 1 to the regulators e and 5 on the upper fermion lines and replace e — > e\ and 8 —* 8\. Sim- 
ilarly for the lower lines we replace e by €2 and 5 by 62 there. Hence, the diagram now depends 
on four variables G(e u e 2 , 81, 8 2 ). Then e 2 G(0, 0, 0) = G(e,0,e,0) + G(e, 0, 0, e) + G(0, e, 0, e, 0) 
+ G(0, e, 0, e) and each of these four diagrams can be reduced to chain integrals which completes 
the evaluation of this graph. 

The calculation of the right hand diagram of figure ^ is more involved. In momentum space 
it corresponds to 

W(A,e)Pj- n (p) 

n(A,e) = ^LUa , (4-H) 

where the transverse projector arises from the fermion loop and the nontrivial piece, H U p, comes 
from the two loop master diagram. We begin with the analysis of the analytic structure of the 
two loop integral. First, it is superficially divergent so one has 

iL V p{A,e) = ri U p—£— + regular terms . 

Second, for the same reason as we discussed for the previous graph n^ p (A,e) — > as e — > 0. 
Third, the residue at the A pole is independent of e. Indeed, since the divergence is logarithmic 
the residue does not depend on the extermal momenta. Then, choosing a new route for the 
external momenta flow, for example along the vertical line, one can see that all the dependence 
on e disappears. Since H vp (A,0) = 0, we conclude that R = and therefore the diagram is 
finite. Moreover, it is not hard to check that the terms proportional to r] up are even in e, and 
therefore we have 



IMA, e) 



e 2 Vup A + P -^(eD + e 2 B + eAC) 



(4.12) 



For e = A the diagram with contracted indices can be reduced, in coordinate space, to the sum 
of chain integrals and one with a unique vertex of the Gross-Neveu type. Next, it can be checked 
that g vp I\. yp {A 1 A) = which results in the following constraints on the coefficients 

D = , 2pA + B + C = . 

We recall that due to the presence of the transverse projector in ( |4.11| ) we are only interested 
in the coefficient A. However, it is easier to calculate coefficients B and C and find A from the 
above constraint. To determine C we repeat the technique used for the evaluation of the previous 
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diagram. In other words we set e — ► ei( 2 ) for the upper and lower lines, respectively. Then, C 
is equal to the coefficient of the eA term in the sum of the two diagrams Yl up (A, e\ = 0, e 2 = e) 
+ Il^p(A, ei = e, £2 = 0), which are easy to integrate. To deduce B we first transform the master 
diagram into coordinate space, giving 

( x 2\m-2A Kvp( x > A, e) • 

It is easy to check that for A = 0, only the term proportional to B survives in II with U up (x, A, e) 
~ B{j] vp — 2fix u x p /x 2 ) + 0(A). Next, we use tetrahedral symmetry whereby we add the new 
line, x u x p /(x 2 ) 1+K , to the diagram to obtain a vacuum graph. The integration over x yields 
the pole in k which is independent of which line we place the regularizing parameter, k, [pSfl . 
Thus we can place the regulator on the vertical line and change the order of integration where 
we integrate over the upper vertex last. After some algebra the new diagram can be reduced to 
chains which allows us to determine the coefficient B, and hence the coefficient A which we are 
interested in. 

Clearly, the evaluation of the three loop contribution to the gluon propagator is a tedious 
exercise. To ensure that we have determined it correctly, aside from the checks we will discuss 
later, we have also undertaken to calculate it by another method. As this is equally as long 
an exercise we will briefly summarize the main steps. It is based on the original method of 
subtractions of [|l2| but differs from the one outlined above in that the original Feynman diagram 
is broken up into a sum of scalar integrals by taking spinor traces and the relevant Lorentz 
projections. Although this results in a large set of integrals the majority of them are in fact 
reducible to simple chain integrals or graphs involving the two loop self energy master diagram. 
The latter can readily be evaluated by uniqueness methods. The remaining three loop graphs 
fall into two classes. Either they are divergent, and therefore their simple pole and finite part 
with respect to A need to be determined, or they are fully A-finite. In the former case we 
were able to apply integration by parts and related methods to again reduce them to integrals 
which were chains, master diagrams or additionally three loop integrals which were evaluated 
by subtraction methods to the finite part. The purely A-finite integrals, in fact, represented 
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Figure 8: Finite three loop diagrams contributing to gluon self-energy. 

the most difficult to determine and were of the form given in figure ||] where n = 0, 1, 2 or 
3 and k, I and q are the loop momenta with p the external momentum. In this notation the 
negative integer exponent corresponds to a numerator factor in the integrand, ((k — l) 2 ) n . The 
n = case had been treated previously in [12]. To determine the other graphs we considered a 



related diagram which involved fermion lines constructed in such a way that when the n traces 
were computed the original scalar diagram of figure || emerged. The additional integrals which 
accompanied this were again either chains or two loop self energy graphs, aside of course from 
the new fermionic graph. This was evaluated by applying a Fourier transform to map it to 



20 



coordinate space before taking the spinor traces. However, to avoid intermediate singularities in 
this step which arise from anti-uniqueness we introduced a temporary analytic regularization, 
5, in the exponents of the propagators. The advantage of transforming to coordinate space is 
that the resulting scalar integrals all have a vertex where the scalar uniqueness relation can be 
applied to leave chain integrals. For completeness we note that the values of the integral of 
figure |8| are 
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for n = 2 and 
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for n = 3. Including these with the other values we have verified that the same value for the 
three loop gluon self energy is obtained in the Landau gauge. 

This completes our detailed discussion of the evaluation of the three loop self energy diagrams 
and we now collect all the contributions and record the final values for all the relevant Feynman 
diagrams. First, we consider the gluon polarization tensor. To fix normalization we write down 
the expressions for the propagators in momentum space 

1 „ GM 2A 
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where the amplitude G in the gluon propagator reads 

G = - . (4.15) 

2r 2 ( M )r(2 - m) K ' 

Further, the contribution from the i-th. graph of figure |] to the polarization tensor can be 
written in momentum space as 
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fJtU 



(4.16) 



Here P^, P" are the transverse and longitudinal projectors and rii is the number of the gluon 
lines in the diagram. The calculation yields for the coefficients B % 
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As can be seen they are all finite and their sum vanishes, Y^,i B l = 0- 

Instead of recording the values for A 1 we list the values for the coefficients II*, which are 
related to the former as 

IT - B i 



A 1 



(2A* -1) ' 

and, as is easy to see, they are given by the trace of the corresponding contributions to the 
polarization tensor, II* = IP,/. We find 
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Here 
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where VK 2 -) = (hir(x))'. Taking into account (|4.16| ) the renormalized gluon propagator is 
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and the functions IIq and IE® are defined as 
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Finally, r/o is twice the anomalous quark dimension in Landau gauge which is recorded in 

0-2)(2 M -l)r(2 M ) 
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(4.28) 



We note that the part of the polarization tensor corresponding to the QED contribution is 
independent of the gauge fixing parameter £ which is consistent with a general analysis. For 



example, see [25]. 



Next we turn to the quark propagator and record the values for the graphs of figure EE First, 
for completeness we write down the value for the graph 




which contributes to the fermion propagator at 0(1/Nf). Straightforward calculations give 
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for the above graph which results in the following anomalous quark dimension at 0(1/Nf) 
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and the renormalized propagator 
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Here and elsewhere we use the notation A = J2k A{k) I 



n 



Below we list the ICR' values of the graphs of figure EJ, which are necessary to determine 
the quark anomalous dimension at Oil/Nj). The value of each graph after application of the 
ICR' operation will be denoted by the respective capital letter. So, for example, applying ICR! 
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to graph (b) of [2] gives ifB. In addition AS 1 ' corresponds to the value of the i-th subgraph of 
figure |||. We find 
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According to ( |3.13| ) the anomalous quark dimension at 0{1/Nj) is given by the sum of the 
residues at simple poles in A multiplied by the number of gluon lines in the diagrams. The final 
result for 772 will be given in section 7. 



5 Schwinger Dyson approach. 

In this section we present the Schwinger Dyson, (SD), formalism to determine the exponent 772- 
There are various motivations for this. One is that in the context of the newer approach of 
the previous sections it will allow an interested reader to compare and contrast the technology 
involved. Though in fact it will not represent a significant amount of extra work given that 
the difficult Feynman integrals have now been evaluated. Second, it provides another check on 
our calculations as it is clear these are quite technical. Third, we would like to discuss some 
distinctions in the treatment of the theories with abelian and nonabelian gauge symmetries in 
the SD approach. Our first comment concerns the issue of choosing the gauge fixing parameter, 
£, entering the gluon propagator. To appreciate the subtlety of this choice we recall that in the 
original construction of Vasil'ev et al for scale invariant theories, such as c/> 4 theory, [|l^], one 
omits at the outset the contribution of the bare propagator in the SD equation, which is less 
singular in the infra-red region compared to the loop corrections. This converts the SD equations 
into self-consistency equations which determine the exponents of the propagators. In respect 
of the discussion of the previous sections concerning QCD and the NATM, scale invariance is 
only present in two gauges. These are the Landau gauge, £ = 0, and that given by £ = £a 
in ( 3.21| ). Therefore, strictly speaking, our subsequent SD analysis will only be valid for these 
particular gauges. By contrast in an abelian gauge theory such as QED the gauge propagator is 
scale invariant in any gauge in the 1/Nf expansion. To understand how the choice of these two 
gauges arises from demanding scale invariance one can repeat the arguments given in section 2. 
Indeed, on general grounds the form of the full gluon propagator in the massless limit will take 
the form 



(5.1) 
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Clearly to have a scale invariant gluon propagator the exponents of the momentum prefactor in 
the transverse and longitudinal pieces both ought to be the same. This can come about in two 
ways subject to the physical restriction of retaining the transverse part. One case is to have no 
longitudinal part which gives rise to the Landau gauge, £ = 0. The other is to match the powers 
of the exponents in each term. This requires the gluon to have zero anomalous dimension, r] 
+ x = 0- Since this exponent is gauge dependent then one can in principle solve this equation 
order by order in large Nf perturbation theory to determine the explicit value of £. It is this 
solution which we denote by £a- It is worth noting that in QED, [p2| , |23|, the photon anomalous 
dimension vanishes in all gauges by virtue of the QED Ward identity and therefore the ansatz 
for the photon propagator was automatically scale invariant. Alternatively another point of 
view can be taken on this problem. One can solve the SD equations in a nonabelian gauge 
theory taking a scale invariant ansatz for the gluon propagator with an arbitrary gauge fixing 
parameter. This parameter is then considered as an input parameter which is to be determined. 
Therefore, it is not hard to understand that solving the SD equations one will find that the 
latter are only consistent for two particular values of £. These will be £ = and £ = ^a- 

In light of these observations we now proceed with the SD calculation but restrict ourselves 
to the Landau gauge. Therefore, the ansatz for the gluon propagator in this calculation is, 



k^k v 



(5.2) 



in momentum space. The gluon propagator exponent is given by 

P = 1 - V ~ X (5-3) 

where x is the quark gluon vertex anomalous dimension. For the quark and ghost fields we 
define their momentum space propagators as 

m = ^ . m _ ^ (5.4) 

where A, B and C are the momentum independent amplitudes. The analogous coordinate space 



amplitudes will be A, B and C respectively. In the SD critical point method, [12], one determines 



the unknown exponent, 77, by representing the two point SD equations at the appropriate order 
which is 0(1/Nj). The representation of the equations form a set of algebraic equations with 
various unknowns one of which is rj whilst the others are combinations of the amplitudes A, B 
and C. Eliminating the latter variables allows one to deduce rj. 




Figure 9: Quark Schwinger Dyson equation to 0(1/Nj). 

To illustrate these points we focus initially on the quark equation as there are additional 
features which need to be considered in the treatment of the gluon equation which will be 
discussed later. If we denote by Si and £2 the values of the respective two and three loop 
graphs of figure |9] in coordinate space then the SD equation at 0(1/N?) can be represented 
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algebraically by 



r(a-l) + S zm 2 (x 2 ) x+A + C F (C F - C A /2)(x 2 ) 2x+2A z 2 Z 1 

+ C F C A (x 2 fx +3A T F N f z 3 i: 2 (5.5) 



where 

09 -A) 
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(5.6) 



X2/J.-20- 1 + 2A) 

corresponds to the first graph of figure || and z = ^4 2 £?. The analytic regularization has been 
included by shifting the vertex anomalous dimension from x to x + A. As in the original method 
of [12] we have chosen the coordinate space representation of the Schwinger Dyson equation. To 
relate to the momentum space version one applies the usual Fourier transform 

1 _ a(a) f e ikx 



{x 2 ) a 2 2a ^ J k (k 2 ) 



\fi-a 



to the version of (5.5) prior to the powers of x 2 being cancelled. The coordinate space propagators 
can be determined from ( |5.2[) and (|5.4j) and are taken to be 



^( x ) = 7~jL > A ^( x ) ~ 
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2/3 x^x u 
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The inverse quark propagator, which is the origin of the first term of (|5.5D , is given by 

x r(q-l)^ 
A(x 2 ) 2 ^- 2a+1 

where 

. . aa(a — u) , . 

r a = -=-r-i . . 5.10 

1 ; tt 2 ^(/x - a)a(a) 



In ( |5.5| ) the quantity m corresponds to the vertex renormalization constant which will absorb 
the divergences arising in the two and three loop corrections. To make this more explicit we 
define their A-dependence as 

^ = § + SJ (5.11) 

where the O(A) terms are not important at this order in 1/Nf. Therefore the quark Dyson 
equation is rendered finite by defining the counterterm formally as 

where we have set 

m - 1 + WE + ■ 2 " t + if + ■ < 5J3 » 

However, this leaves terms in the SD equation involving lnx 2 which would otherwise spoil the 
simple scaling behaviour at the fixed point. To avoid this one defines the vertex anomalous 
dimension to be 
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With these definitions the finite quark Dyson equation to 0(1/ N 2 ) is 
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where the first term contains r)\ and 7/2 and we have expanded Sq as 
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Therefore, from ( |5.15| ) we can determine a relation for 772 in terms of z<i and the finite correction 
graphs, as 



z 2 m rfi 
Z\ 2\x 
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(5.17) 



where r\ = Y^=i Vi/Nf and 771 = 7^ > /(2T F ) in earlier notation. 

For the gluon SD equation there is the added complication of dealing with the transverse and 
longitudinal components of the equation. Whilst this has been dealt with for QED in p3| we 
will recall the important features. In coordinate space the inverse gluon propagator is, fl32|, 
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We recall, |3^, |23|], that this is deduced by inverting the gauge dependent momentum space gluon 
propagator on the physical transverse subspace before Fourier transforming back to coordinate 
space. Therefore, we can formally represent the coordinate gluon SD equation as 
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(5.20) 



where 111 and Hi relate to the coordinate space value of the two loop diagram of figure pi] 
and II2 and H2 correspond to that of the three loop one. The final term of ( 5.20|) is the ghost 
contribution from the ghost graph of figure [[(]. However, its treatment requires special care. 
Although we have constructed the gluon equation in coordinate space, the ghost contribution is 
determined from its evaluation in momentum space since the ghost couples via a derivative to 
the gluon. This explains the appearance of the amplitude combination y = C 2 B in fl5.20p . In 
particular the ghost graph of figure |l0| is 
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Figure 10: Gluon Schwinger Dyson equation to 0(1/ Nf). 



Although the ghost was treated in large Nt in the Landau gauge in |15| where its anomalous 
dimension and vertex anomalous dimension were computed at 0(1/Nf), the value of yi was not 
recorded. To determine it we consider the ghost SD equation of figure [Tl]. In momentum space 
it becomes 

u(/j,--y -1,/j,- /3,/3 + j) 



whence 



= 1 + (2/x-l)yCU 
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where r\ c 
we have 



(2fi - 1)C A 

l=1 rjcj/Nj. Alternatively using the Slavnov- Taylor identity in the Landau gauge 

VI = 7 r( ^ +1)?? \ . (5.24) 
2( A t-2)(2/ i -l)T F V ; 

We note that in ( 5.20| ) the ^-dependence of the final term involves \c- This emerges after using 
the Slavnov- Taylor identity in rationalising powers of x 2 in the initial representation of the SD 
equation. 

Whilst there appears to be two components of (|5.20|) to analyse it turns out that only one 
is important. This is the transverse component when expressed in momentum space using the 
Fourier transform. Although this may appear to neglect information contained in the longitu- 
dinal part of the equation with respect to momentum space it turns out in fact that the sum 
of contributions to the longitudinal projection of ( |5.2C ) in momentum space from the various 
graphs in arbitrary gauge is actually zero. The cancellation in the Cp sector was discussed in 
[23]. For the Ca sector the contributions from the longitudinal sectors of the last three graphs 



of figure 10 sum to zero which can be verified from the explicit values of the graphs. This is also 
true for non-zero £. 



Therefore, returning to ( 5.20| ) and restoring the common factor of l/(x 2 ) 2fM - f3+A in each 



term we transform the equation to momentum and retain only the transverse part. Inverting this 
component to coordinate space yields the relevant part of the gluon SD equation for determining 



= 



Figure 11: Ghost Schwinger Dyson equation at 0(1/Nf). 
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the critical exponents, 
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which now requires renormalization. This is performed in the same way that (5.5) was rendered 
finite. Repeating the procedure here and ensuring that there is no anomalous scaling behaviour 
we find, 
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Although it may seem that ( |5.12| ) and (5.14) will give different values from ( |5.28 ), when the 
explicit values for the residues of the simple poles in A are substituted, they are equivalent. 
Moreover, the value of xi wm agree with ( 5.14j ). Having made the gluon Schwinger Dyson 
equation finite, it is elementary to expand ( 5.26| ) to 0(1/N 2 ) to obtain an expression for z 2 . We 
find 
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Therefore with this value for z 2 we can establish a formal expression for r] 2 in arbitrary gauge 
which depends only on the values of the various graphs of figures [9l |l0| and [ll| by eliminating 
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Z2 in ( 5.17 ). For completeness, we quote the coordinate space values of the integrals which will 
determine 772 from these equations. First, we recall the graphs of the QED sector are 



Hi 



/x(2//-3) 2 r 2 (//) 



"1 
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2(2//-l) 2 (//-2) 



+ 4(/x - 1 
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A 
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(m-i) 

For the three loop graphs we have 

8(2//-l) 2 (//-2) r(2//-l)(//-2)(//-l) 



M(2/x - 3) 



(5.30) 



(5.31) 
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(5.32) 



The longitudinal component is related to 1I2 by 
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-2 



2Il2 + 



(2//-3) 2 // 2 r4(/,)r ?( 



'o 



(2// -!)(// -2) + 



//(//- l)V(l) 



2(2// - 3)o(3 - //) 



(5.33) 



(5.34) 



Therefore, substituting the values of the integrals in the formal expressions for xi an d mi, we 
have explicitly 



mi 
Xi 



C F rjo (2// - 1)(// - 3)CU?7o 



2T F 
CfVo 



8(2//-l)(//-2)Tp 
C^r/o 



2(/x - 2)7> 



(5.35) 
(5.36) 



The relation between the coordinate space variables and the momentum space ones of earlier 
sections is determined from the Fourier transform of the asymptotic scaling forms. For instance, 
it is easy to deduce that the respective quark amplitudes A and A are related by 



A 



ia(a — 1)A 



(5.37) 



(a-1) 

whilst the respective gluon amplitudes are determined from comparing the coefficients of the 
77^2, component. We find 

2(n-0)a(0)B 
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(2// -2/3-1) 



(5.38) 
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Thus, 



which implies 



and 



z = — 



{[i- (3)a 2 {a-l)a{l3)z 
(a -1)2(2^ -2/3-1) 

2zi 



(2/x - 3)I» 



Z2 



Z2 



(2/i - 3) 



(5.39) 
(5.40) 
(5.41) 



(2fi - 3)I» 

upon expanding in powers of l/iVf. Moreover, variables involving the ghost amplitudes are 
related by y = a 2 {~f)a{0)y. Similarly, we can relate the values of the three loop integrals 
expressed in either coordinate or momentum space. For instance, if the value of the three loop 
gluonic graph in coordinate space is 



1 



and in momentum space 



(• x 2)2^-l-2A 
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n „ XfiX u 
+ -2 



fr . ,=, PfiPv 
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Likewise, 
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(5.42) 
(5.43) 



16(/x- l) 5 (2 / u-3 + 2A) 3 a 3 ( / u-l + A)a( / u-3A) ' 
So, for example, the finite parts of the momentum space integrals are given by 
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(2/x-l)(2/x-3)(/x-l)( M -2)/x 2 
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(5.46) 



and 



(2/i-3)( M -2)/i 1 
2(/x-l)(2/,-l) 2 (/x-2) 2 



4(2/x-3)(/x-2)- 



0(M) 



a(3 — jj) 

(4// 5 - 6// 4 - 13/x 3 + 40/r 2 - 40/x + 10) 
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(5.47) 



(5.48) 



31 



6 Mass anomalous dimension. 



Having presented the formalism to calculate the quark anomalous dimension at 0(1/Nf) we are 
now turn to the quark mass anomalous dimension at the same order. Perturbatively the mass 
anomalous dimension is known to four loop accuracy. Thus comparison of results obtained in 
both perturbation theory and the 1/iVf expansion will provide a non-trivial test on the validity 
of our results since, for instance, only three loop diagrams are considered at 0(1 /Nf) in the 
large Nf formalism. To complete the calculation we discuss some observations here which allow 
us to relate the values of most of the graphs contributing to the mass dimensions at 0(1 /Nf) 
to those used for the quark anomalous dimensions. It transpires that this greatly reduces the 
amount of calculation and will provide a method for tackling the evaluation of the anomalous 
dimensions of other operators. 

As is well known the mass anomalous dimension coincides with the anomalous dimension of 
the mass operator ipip, 

7m = Ifrfi + . (6.1) 

Here 7^ is the quark anomalous dimension and the RG function 7^ is expressed via the renor- 
malization constant of the mass operator by 

7^ = Md M lnZ^ , [4>^] R = Z^Z^ 2 [VWo]o • (6.2) 

To determine the renormalization constant Z^ at second order in the 1/Nf expansion one 
has to compute the divergences of the one-particle irreducible 2-point Green function with the 
insertion of the mass operator. The corresponding diagrams can be obtained from the diagrams 
contributing to the quark self-energy of figure ^ by the insertion of a Tpij) vertex in the quark 
lines connecting the external vertices. At the level of diagrams this equivalent to the insertion 
of the unit operator in a fermion line 

As was discussed in section 4 the differentiation of the propagator diagrams with respect to 
Pa, where p is the external momentum, results in the same set of the diagrams as for the 
mass operator. The only difference is that the inserted vertex in this case is (—7^). We shall 
exploit this fact and show below that provided the Landau gauge is used then the corresponding 
contributions to Z\ = Z^ and Z^w, from all diagrams from both sets, except for two which arise 
from graphs (b) and (c) of figure |2] when the insertion of the new vertex is in the middle quark 
line, are related to one other by 

Z f = . (6 .3) 

So given that 7^, is known then to determine the mass dimension at 0(1 /Nf) it is sufficient 
to compute the contributions to Z\ and Zj^ from the two extra diagrams. However, their 
calculation does not lead to any difficulties using the methods discussed in section 4. 



We now proceed to the proof of the relation ( |6.3p . First, we use the freedom we have to 
choose the external momentum flow in a diagram arbitrarily when calculating the contribution 
to the renormalization constants by directing the external momentum flow through the inserted 
vertex and out through the nearest external vertex. After this all the diagrams in question take 
the form shown in figure 12. All that one needs to know about the coloured block in figure ^ 
is its Lorentz structure, 

r,(* A) = (j v A(A) + ^B(A)) , (6.4) 
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Figure 12: External momenta routing in the quark 2-point function with an operator insertion. 
The large black dot stands for the Vi = {I, (— 7 M )}. 



where rii is the number of the gluon lines in the block and the exact form of the functions 
A{A)andB{A) are irrelevant for our present discussion. Since the gluon propagator which joins 
this generalized vertex is purely transverse in the Landau gauge the second term in fl6.4| ) will 
not contribute to the final answer. Thus the momentum integral for the diagram in figure ^ 
can be cast in the form 



A(A) 



d 2 ^q l p {i>-i)Vil v jP^{q) 
(2vr) 2 ^ (p - q) 2 ( g 2) Al +(n i +l)A 



(6.5) 



We are interested in the pole terms of this integral only, and bearing in mind that after the 
subtraction of the divergent subgraphs, which also have the same form, the latter is independent 
of p. Therefore one can rewrite (|6.5| ) as 



A(A) Vi 



1 

2~M 



(2vr) 2 M (p - q) 2 ( g 2)A»-M»i+l)A 



(6.6) 



where we have used the fact that up to p-dependent terms one can replace q p q u by q 2 /(2p) 
in the integrand. It is easy to deduce that the combination in the prefactor in ( |6.6[) is equal 
to — (2/u — 1) for V = I, and to (2/i, — l)(/i — 2)/// for V = — 7 M . Then one can immediately 
conclude that irrespective of the exact form of the constant ^4(A) which corresponds to the 



structure of the coloured block in figure 12, the contributions of the corresponding diagram to 



the renormalization constants Z\ and are indeed related by equation (16^) . To complete the 
calculation of 7 m one needs to compute the contributions from the graphs originating from those 
b, c in figure These, as well as the determination of the graphs when the gluon propagator is 
taken in an arbitrary gauge, can readily be performed. 



7 Discussion. 



We are now in a position to assemble our results. First, using the formalism of section 3 we 
record that the (i-dimensional expression for the quark anomalous dimension is 
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£(/i 2 -4^ + 2) £ 2 (2 M 3 - 11/i 2 + 12// - 4) 



2( M -1)(// 



4(a*-1)(m-2)(2m-1) 



(7.1) 



We note that as an initial check on the result setting Cp = 1 and Ca = one obtains the 
expression for the quark anomalous dimension of the abelian Thirring model. Again, in agree- 
ment with general results, [25], we find that in this case 7! is independent of the gauge fixing 
parameter. Further, collecting all the results for the quark mass anomalous dimension, we find 
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/x 2 (2 M 
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4(/i-2)( M -l) 



-l)(2 M -3)(/i-2) 

C F C A r/ 2 
( M - 2) 2 (2 M -1) ' 



(7.2) 



where J m (9*) = Y^°=ilm,r/ri r ■ We have performed the calculations with a non-zero gauge 
parameter to check that it cancels exactly since the mass dimension is a gauge independent 
quantity. Aside from the internal checks on the integration which we discussed earlier, there are 
other checks on the correctness of (7^1) and (|7.2j). First, the terms of ( |7.1[) and (|7.2j ) involving 
Cp have both been derived previously in QED in |23| and [24] respectively. However, a more 



stringent check rests in the relation of the results with the critical renormalization group equation 
since the RG functions evaluated at ought to agree with (|7. 1|) and (|7.2j) at the same order in 



1/Nf and 0(e ). It is elementary to verify that our results are in agreement with the Landau 
gauge anomalous dimension of [||, [35f| . Also the result for j m (g*) is consistent with || [ll], ^, 

Having verified this agreement we can now determine new information on the structure of 
the renormalization group functions in principle to all orders in the strong coupling constant. 
However, we will record values of the coefficients to six loops by writing 

25 
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(7.3) 



for the Landau gauge and 
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(7.4) 



for the mass dimension where the order symbol means that contributions which are third order 
in 1/Nf and 0(e 7 ) at criticality are ignored. Although the full four loop MS /^-function is only 



34 



available it may seem that it is not possible to decode the Landau gauge information in (|7.1[) 
and ( ff^ ) beyond four loops since this requires a* to 0(e 6 ). However, recalling that we are 
computing within the 1/Nt expansion we can use the information contained in the QCD (5- 
function exponent u = — (3'{a*)/2 which has been determined at 0(1/Nf) in ]l6|]. From 
we deduce that the relevant corrections are 
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Therefore, with this value we deduce the new exact contributions to ( |7.3| ) and ([/ 
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(7.8) 



We have given the explicit expression for C41 in order to compare with the full MS gauge cal- 
culation of which was carried out for the particular colour group SU(N C ) and note exact 
agreement. Further, we note that in [19] it was pointed out how the asymptotic Pade approxi- 
mant predictions of 36] for rn-50 and m^i compared with the exact values we have determined. 
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Although we have concentrated on the relation of our results to four dimensional perturbation 
theory, we can also quote values of the exponents in three dimensions. These will be important for 
alternative critical point investigations of the non-abelian Thirring model and QCD equivalence, 



such as the lattice. For the wave function, we find from (7.1) 



4Cf(3£-2) 
i] = „ 0m — — + 



8C F 



9k a T f NJ 
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However, for the physically more interesting mass exponent we deduce, 
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(7.9) 



(7.10) 



It would be interesting to see how the numerical values of these exponents compare with the 
lattice. Indeed the domi: 
arises from the C A term. 



lattice. Indeed the dominant contribution from the 0(1/N 2 ) correction in the latter exponent 



In conclusion we have provided an elegant alternative to computing perturbative information 
on QCD at a new order in the large Nt expansion. Although the evaluation of the Feynman 
integrals has formed a significant part of the discussion, is technically quite involved and provides 
useful techniques for future massless integral evaluation, we believe we have demonstrated that 
the computation of other renormalization group functions in QCD is viable in principle. An 
example of this is the twist-2 operators which arise in deep inelastic scattering. Furthermore, 
the relation of the NATM to QCD has been put on a firmer ground. Indeed we believe that this 
critical relation between both models deserves further investigation. 
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A Basic integration rules. 

In this appendix we give formulas for the Fourier transformation, integration of chains and 
uniqueness relations for the vertices of different type. The following notation is used. We define 

n 

a(xi,...,x n ) = Y[a(xi) 
i=l 

where o(a) = r(a')/r(o!) with a' = (fi — a). The Fourier transforms of various functions are 
given by 
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The integration of the basic chain with various tensor generalizations are given by 
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(A.4) 



Next we list the uniqueness relations for the three type of vertices which arose in our calculations 
where the appropriate uniqueness value for each vertex is recorded in parentheses beside each 
rule. 

Boson vertex (a + j3 + 7 = 2p). 
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Gross-Neveu vertex (a + /3 + 7 = 2/i — 1). 




Vector- fermion vertex (a + /3 + 7 = 2/i — 1). 




Here 



(2 



[z — x) 2 (y — x) 2 

As usual the line with an arrow with index a denotes the fermion propagator, $/(x 2 ) a , a simple 
line corresponds to the boson propagator l/(x 2 ) a and a wavy one is used to denote the conformal 
propagator (4.4). 



B Renormalization constants. 



In this appendix we collect the values of the renormalization constants Z^ and the basic RG 
functions at 0(1/NA. The calculations are rather straightforward and so we only list the results. 
We found 
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and (7 is defined in ( |3.12| ). Using d3.13[ ), ( |3.14| ) and ( |Q| ) one finds 
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For /3a we obtain 
/3a = A( 7 A + 27c-74) 
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One can see that A = 1 is indeed a zero of the beta function. Also, for £4 we find 
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and /3a is simplified in this gauge to 
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Finally, for completeness we record the value of the amplitude combination z required in the 
SD formalism at 0{l/Nf). It is 
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9(m) 



r(/i + l)r?o 
2(2^-1)^-2) 



(B.9) 



+ 



4(2^-1)2(^-2) 



30(/i) + 



(8^ 5 - 92// 4 + 270/i 3 - 301/x 2 + 124/x - 12)*(/x) 



- tf 2 (/x) - $( M ) + 



+ 



2fi(fi - l)(2/x - l)(2/z - 3)(/x - 2) 
(p 2 + 2 M -2)£ 



4(/i-l) 2 2(/i-l)V 
(16/x 7 - 120^ 6 + 420 / t/ 5 - 776/x 4 + 742fi 3 - 349^ 2 + 84/x - 12) 
2(2// - 1)(2// - 3)(/x - l) 2 (/i - 2)/i 2 

where the gluon propagator in the SD formalism is taken as 

B 



(B.10) 



h h 



(B.ll) 



in momentum space. These will be useful in the calculation of gauge independent critical expo- 
nents at 0{1/Nj). For example, if one chooses a particular gauge to simplify those calculations, 
such as the Feynman gauge used in standard perturbative calculations, then one will require the 
value of the relevant variables for a general gauge parameter. 
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